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Abstract 

A dynamic phase transition (DPT) with respect to the period P of an applied alternating mag- 
netic field has been observed previously in numerical simulations of magnetic systems. However, 
experimental evidence for this DPT has thus far been limited to qualitative observations of hystere- 
sis loop collapse in studies of hysteresis loop area scaling. Here, we present significantly stronger 
evidence for the experimental observation of this DPT, in a [Co(4 A)/Pt(7 A)]3-multilayer sys- 
tem with strong perpendicular anisotropy. We applied an out-of-plane, time-varying (sawtooth) 
field to the [Co/Pt]3 multilayer, in the presence of a small additional constant field, H^. We then 
measured the resulting out-of-plane magnetization time series to produce nonequilibrium phase 
diagrams (NEPDs) of the cycle-averaged magnetization, Q, and its variance, a 2 (Q), as functions 
of P and Hf,. The experimental NEPDs are found to strongly resemble those calculated from sim- 
ulations of a kinetic Ising model under analagous conditions. The similarity of the experimental 
and simulated NEPDs, in particular the presence of a localized peak in the variance a 2 (Q) in the 
experimental results, constitutes strong evidence for the presence of this DPT in our magnetic 
multilayer samples. Technical challenges related to the hysteretic nature and response time of the 
electromagnet used to generate the time-varying applied field precluded us from extracting mean- 
ingful critical scaling exponents from the current data. However, based on our results, we propose 
refinements to the experimental procedure which could potentially enable the determination of 
critical exponents in the future. 

PACS numbers: 64.60.Ht, 75.70.Cn, 75.60.Ej, 75.10.Hk 
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I. INTRODUCTION 



A dynamic phase transition (DPT), in which the dynamical behavior of a nonequilibrium 
system changes in a singular way at a critical value of a system parameter, can provide 
insight into the often complex behavior of such systems. DPTs of various kinds have been 
identified and studied in chemical,- charge- density wave,- and superconducting systems.- 
The DPT of interest here was first identified in computer simulations of magnetic systems.- 
When a bistable magnetic system in its ferromagnetic phase is exposed to an alternating 
magnetic field, the response depends strongly on the period of the applied field. At low 
values of the field period, the system effectively cannot respond to the rapidly changing 
field, and its magnetization oscillates in a restricted range around one of its two nonzero 
(magnetized) values. At high values of the period, the magnetization can follow the field, 
resulting in a square, symmetric hysteresis loop. This behavior suggests the cycle-averaged 
magnetization Q (with a value near ±1 at low field period, and near at high field period) 
as a candidate for a 'dynamic' order parameter. It has indeed been shown computationally 5 
that at a critical period P c , there exists a second-order phase transition with respect to Q, 
with critical exponents consistent with those of the equilibrium Ising transition. For further 
description of this DPT and the numerical evidence for the second-order phase transition, 



see Refs. 
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Also associated with the appearance of a non-zero value of Q below the critical period P c 
is a significant decrease in the hystersis-loop area. Since first being identified, this DPT has 
received a great deal of attention in numerical simulations^ 1 ^ and analytical work.— 
However, to date there have been only tentative experimental indications of its presence, 
principally the collapse of the hystersis loop area with decreasing field period, in several 
studies of hysteresis scaling properties in magnetic thin films.— 'i 5 - 

It is important to understand the extent to which this DPT, which has been well stud- 
ied computationally, can be realized in an experimental system. To achieve this, one must 
investigate the quantitative behavior of the dynamic order parameter itself. To this end, we 
identified an experimental system, an ultra-thin [Co(4 A)/Pt(7 A) ^-multilayer, which pos- 
sesses strong perpedicular anisotropy and which exhibits substantial similarities 16 to a two- 
dimensional kinetic Ising model,— ii£ii2i22i2ii22i22i2i for which the DPT is well-documented.-! 6 - 
We present here the first quant it it at ive data on the behavior of the dynamic order param- 
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eter in this Ising-like experimental system, and we compare it to the behavior observed 
in kinetic Monte Carlo (MC) simulations. Recent computational work— by our group on 
the two-dimensional kinetic Ising model identified the cycle-averaged magnetic field, Hf,, 
as a significant component of the field conjugate to Q, and established the existence of a 
fluctuation-dissipation relation (FDR) between the susceptibility dQ/dH^ and the variance 
<J 2 (Q) in the vicinity of the critical period. Applying those results here, we use the variance 
cr 2 (Q) as a proxy for the susceptibility dQ/dHt in our comparisons of experimental and 
computational results. 

In order to establish the presence of a phase transition unequivocally, it is necessary to 
demonstrate power-law scaling with well-defined critical exponents. While this often can 
be done straightforwardly enough in computational work, experimentally it requires precise 
measurements at carefully controlled values of the relevant thermodynamic variables. This 
was achieved for the equilibrium Curie transition in an experimental thin film only fairly 
recently— In the present case of a [Co/Pt]3 multilayer driven by an alternating magnetic 
field, the nonlinear, hysteretic nature of the electromagnet used to generate the alternating 
field resulted in small fluctuations of the cycle- averaged field, Hb, rendering it impossible to 
study the experimental system at if& = 0, as done in all but the most recent^ computational 
studies. Moreover, it was not possible to establish scaling relations with statistical signifi- 
cance from our data, despite achieving an applied field and bias field amplitude precision of 
0.5% and 0.1%, respectively, in our measurements. Nonetheless, our data for Q and <J 2 (Q) 
are consistent with power- law scaling near the critical point. 

In spite of this, we argue that the similarity of the behavior of the cycle-averaged mag- 
netization Q in the multilayer to its behavior in the kinetic Ising model, coupled with an 
explanation of the differences which do exist in terms of known physical properties of the 
multilayer, provides strong evidence for the presence of this DPT in the multilayer. In 
addition, we propose refinements in the experimental technique which may produce in the 
future the precision required to investigate power-law scaling. However, peaks in response 
functions also consitute evidence for phase transitions, and our observation of a peak in the 
variance cr 2 (Q) represents significantly stronger evidence for this DPT in an experimental 
system than that previously published.— 1 ^ 

The importance of the present study is twofold. First, it provides new insight into the 
dynamic behavior of an ultrathin magnetic film system whose properties are important 
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to the development of future generations of ultrahigh density information-recording tech- 
nologies. Second, it provides valuable experimental input to theorists' efforts to develop a 
comprehensive understanding of phase transitions in far-from-equililibrium systems. 

This paper is organized as follows. In Section [III we describe the experimental setup 
and procedure. In Section IHH we present our experimental results, consisting of directly 
measured hysteresis loops and magnetization time series, as well as nonequilibrium phase 
diagrams (NEPDs) used to characterize the behavior of the dynamic order parameter and 
its fluctuations. In Section HVl we compare the experimental results to those of computer 
simulations of a kinetic Ising model and argue that the similarities provide strong evidence for 
the presence of the DPT in the multilayer system. In Section |V] we present our conclusions, 
as well as suggestions for further experimental work. 

II. EXPERIMENTAL SETUP 

The ultra-thin [Co (4 A)/Pt(7 A)]3-multilayer samples were prepared by low-pressure 
(3 mtorr Ar) magnetron sputtering onto ambient temperature silicon-nitride coated Si sub- 
strates. We deposited 200 A Pt as a seed layer, and the samples were coated with 20 A Pt 
to avoid contamination. These ferromagnetic multilayers have an easy axis along the surface 
normal and strong uniaxial anisotropy, resulting in high remanent magnetization and square 
out-of-plane hysteresis loops.^ 7 X-ray diffraction confirmed a (111) crystalline texture, with 
a lateral crystallographic coherence length from several tens to several hundreds of nm. 

Due to the strong ferromagnetic interlayer coupling, mediated by the Pt between ad- 
jacent Co layers,— the entire multilayer acts as a single magnetic film with strong uniax- 
ial anisotropy. While out-of-plane magnetized ferromagnetic films tend to form equilib- 
rium domain structures to reduce magnetostatic interactions, in ultra-thin films the en- 
ergy reduction resulting from domain formation is extremely small, so that this effect is 
strongly suppressed.— \y e measured the effective anisotropy field of the multilayer to 
be Hk ~ 6 kOe, more than an order of magnitude larger than the fields used in our experi- 
ment. The spins should therefore remain strongly collinear perpendicular to the film. These 
experimental facts, combined with theoretical and experimental 24 evidence that ultrathin 
films with uniaxial anisotropy are in the same universality class as the equilibrium Ising 
model, suggest that the two-dimensional kinetic Ising model can serve as a useful model for 
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the nonequilibrium behavior of our multilayer. 

To study magnetization reversal in the multilayer, we measured the polar magneto-optical 
Kerr effect (MOKE), which is proportional to the magnetization along the surface normal. A 
time-dependent magnetic field along the surface normal was provided by an electromagnet. 
In addition, for reasons discussed in Sec. IHH a small constant additional magnetic field - 
the 'bias field', Hb - was applied by separate means. The total (time- varying) field actually 
experienced by the multilayer was monitored by a Hall probe. Data were recorded for two 
multilayer samples, A and B. The electromagnet was driven by a computer-controlled bi- 
polar power supply to provide a saw-tooth field sequence for several values of the period 
between P = 8.7 s and 62.3 s (sample B: between P = 7.6 s and 39.6 s). For each sequence, 
we first measured two complete, saturated hysteresis loops using a large field amplitude 
H s = 740 ± 5 Oe (for both samples A and B). Next, in the actual measurement sequence, 
the amplitude was lowered to H = 366 Oe (sample B: H = 344 Oe), and 49 identical cycles 
were run. By monitoring the field sequence, we determined that the variation of Hq during 
the 49 cycles was less than 1%. Finally, the amplitude was increased back to H s to trace 
out another complete hysteresis loop. 

The two previous experiments^ 1 ^ which observed hysteresis loop collapse with decreasing 
field period were conducted on 3 monolayer (3 ML) Co/Cu(001) ultrathin films with in- 
plane uniaxial anisotropy^ and on 2-3 ML Fe/W(110) ultrathin films with perpendicular 
uniaxial anisotropy.~ These experiments also employed an electromagnet to generate the 
time-varying magnetic field, a Hall probe to record the magnetic field at the location of the 
sample, and the MOKE effect to measure the sample magnetization. In our experiment, 
however, two additional steps were taken to enable a more thorough investigation of the 
experimental behavior of the cycle-averaged magnetization, 



previously identified as the order parameter for the DPT in the kinetic Ising model.- (The 
index % in Eq. (JTJ) denotes the number of the field cycle.) 

First, as mentioned above, a significant component of the field conjugate to the dynamic 
order parameter was identified in a recent computational study^ as the cycle-averaged value 
of the magnetic field, or 'bias field', 




(1) 




(2) 
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In Ref. [25|, with Hb taken as a small constant field superposed on a time-varying square- 
wave field, the conjugate field scaling exponent 5d in the DPT was shown to agree with the 
equilibrium Ising field scaling exponent 5 e = 15 to within computational error. The bias field 
was found to have a signficant effect on the dynamic order parameter Q, especially near the 
critical period where the susceptibility dQ/dH^ has its maximum. In a computer simulation, 
one can study the DPT in precisely zero bias field without difficulty In experiment, however, 
due to nonlinearities in the electromagnet used to generate the time-varying field, the actual 
bias field Hf, experienced by the sample will inevitably exhibit finite values, which may even 
fluctuate slightly during the measurement sequence. We therefore chose to run experiments 
at a series of different values of the bias field, which were measured and controlled carefully 
using the Hall probe, 26 at each value of the field period studied. 

Second, a recent study of thicker [Co/Ptjso multilayers 2 ^ found that hard-to-reverse, resid- 
ual bubble domains could persist beyond an apparent saturation field, and could serve as 
nucleation sites during magnetization reversal when the field was subsequentely reversed. 
In the experiments reported here, the two complete loops at high saturating field before the 
measurement sequence ensured that the samples began each run in a well defined magnetic 
state, and that the effects of residual bubble domains at the lower field magnitude of the 
measurement cycle were consistent across different data runs. In addition, comparison of 
the complete loops before and after the measurement sequence served to verify the stability 
of our experimental setup during each run. 



III. EXPERIMENTAL RESULTS 



The hysteresis loops and magnetization time series in samples A and B had very similar 
characteristics, so we show only those of sample A (the NEPDs for both samples A and 
B will be presented, however). Figure [T] shows the normalized magnetization, m, vs H for 
periods P = 16.2 s and P = 38.1 s in sample A. The initial and final complete loops exhibit 
full remanent magnetization in a single- domain state, and a sharp reversal region. The loops 
in the measurement sequence, at field magnitude H , are also square with a sharp reversal 
region. However, a sharp suppression of the magnetization reversal process is visible below 
a pinning field H p = 288 ±2 Oe (sample B: H p = 260 ± 6 Oe).— The pinning phenomenon 
has been studied quantitatively in multilayer films, and it is understood to result from a 
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distribution of energy barriers to domain wall motion resulting from lattice defects and 
variations in the film thickness.— vM£L£li21 

The Kerr microscope image in Fig. [2] illustrates the domain reversal pattern in a sister 
sample, at magnetization m ~ 0.9 after nucleation from a positive saturated state. It shows 
clearly that, above the pinning field H p , the reversal within our 1.0 mm 2 MOKE laser spot 
occurs by a multidroplet (MD) process^ of nucleation and growth of approximately circular 
domains. This reversal pattern is similar to that previously observed in ultrathin, single- 
layer CoPt films,- 5 and can be contrasted with the reversal pattern in thicker [Co/Ptjso 
multilayers,— in which larger magnetostatic interactions lead to the formation of stripe 
domains. The MD process is quite similar to that occurring near the DPT in the kinetic 
Ising model^ indicating that, if the effect of the pinning field is properly accounted for in 
the analysis, the kinetic Ising model can serve as a reasonable model of the multilayer. 

As described in Section|TTl given the important role of the bias field Hb as shown in Ref.l25. 
we have studied the behavior (the average value and the variance) of Q in the multilayer 
with respect to both P and Hb, and compared this to the corresponding behavior of Q 
in the two-dimensional kinetic Ising model. The two previous experimental reports of the 



collapse of the hysteresis loop with decreasing field period concentrated principal 
decrease of the hysteresis loop area (cf. Fig. 1 of Ref. [14|, and Figs. 1 and 2 of Ref. 



y on the 



15|). The 



cycle-averaged magnetization was also plotted, in Fig. 3 of Ref. [14], where its magnitude was 
shown to become nonzero and then to increase as the applied field amplitude is decreased at 
a constant value of the field period. However, the bias field, which has a strong effect on the 
cycle-averaged magnetization near the dynamic phase transition in computational studies^ 
was not carefully controlled, and the variance in the cycle-averaged magnetization was not 
recorded. 

From an initial value of Hb = —3.9 Oe (sample B: —7.4 Oe), the bias field, as generated 
separately from the time- varying field of the electromagnet, was increased in steps of 0.5 Oe 
(sample B: steps of 1.1 Oe). The actual value of the bias field during the 49 measurement 
cyles was determined, using data from the Hall probe, to fluctuate at the 0.1% level of the 
applied field amplitude, H Q . Specifically, the actual bias field could be characterized by a 
confidence interval (H b - AH b /2, H b + AH b /2), with AH b = 1.0 Oe. 

In Fig. E(a), magnetization time series for sample A at P = 16.2 s are shown in a strong 
negative, weak negative, and weak positive bias field. The corresponding time series, Q$, of 
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the cycle-averaged magnetization are plotted in Fig. [21(b). For both negative Hb values, the 
cycle-averaged value of the magnetization (and in fact of any quantity) settles eventually 
into a non-equilibrium steady state (NESS), so that Qi is observed in Fig. 121(b) to fluctuate 
around an average value (Q) < 0. (The fluctuations of Qi in the NESS arise from variations 
in the minor loop amplitude, as seen in Figfjja), which are caused principally by thermally- 
induced fluctuations in the extent and timing of nucleation events.) The transition to the 
NESS with (Q) < takes longer in the weak negative Hb than in the strong negative Hb, 
and it does not occur at all in the weak positive H b . At the longer period P = 38.1 s, in 
contrast, the time series adjust quickly to a NESS with a small (Q) value of the same sign as 
Hf,, as shown in Figs. [21(c) and [21(d). This behavior will be compared to that in the kinetic 
Ising model in the next section. 

To characterize the behavior of the cycle-averaged magnetization in the multilayer more 
fully, we measured its average (Q) over the field cycles for which the system was in a NESS, 
for a range of values of Hb at each period P. We will refer to the resulting contour plot 
of (Q) vs P and Hb, shown in Fig. IU(a), as a non-equilibrium phase diagram (NEPD), so 
named in analogy to plots of thermodynamic relationships in equilibrium phase diagrams. As 
described in Section [J the recent confirmation^ of a FDR near the critical point of the DPT 
in the two-dimensional kinetic Ising model justifies the use of the variance as a proxy for the 
susceptibility in evaluating evidence for the DPT. The NEPD of the variance cr 2 (Q) within 
the NESS is shown in Fig. H](b). The corresponding NEPDs for sample B are presented in 
Fig. [5j To determine in which field cycle the system should be considered to have entered a 
NESS, we used the following procedure. The first 10 field cycles were discarded to minimize 
initial transient effects. We then identified the first subsequent field cycle, n, for which the 
slope of the remaining values {Qi, i > n} could not be statistically distinguished from zero. 
If more than 15 values of Qi remained, the mean (Q) and variance <J 2 (Q) of these values 
were calculated; otherwise, no data were recorded on the NEPDs for that pair of values 
(P, Hb). This procedure ensured that, for example for the measurement sequences in strong 
negative and weak negative bias fields in Fig. [3](a), the field cycles making up a transition 
to the NESS with (Q) > were not included as part of the respective NESS's. 

In Fig. [U(a), the boundary between regions with negative (Q) and those with positive 
(Q) occurs, for the lower periods, at a slightly negative bias field, Hb ~ — 1 Oe. Simi- 
larly, the maximum of cr 2 (Q) occurs, for the lower periods, at a slightly negative bias field, 
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Hf, »s —1.5 Oe. This appears to contradict the expected M — H inversion symmetry of a 
ferromagnetic system, which predicts that the NEPDs should be symmetric with respect to 
the bias field Hb- However, it can be explained by the presence of residual, hard-to-reverse 
bubble domains, with coercivities H c in the range H < \H C \ < H s , as observed recently 
in thicker [Co/Ptjso multilayers. 27 Since the final saturating field before our measurement 
sequence begins is positive, such residual domains (if present) would remain positively mag- 
netized during the measurement sequence, in which the field magnitudes < Hq would 
be too weak to reverse them. The bubble domains would then serve as nucleation centers for 
reversal on the increasing branch of the hysteresis loops during the measurement sequence, 
lowering the nucleation field on this branch. This effect can indeed be observed in Fig. UK a), 
in which the nucleation field H n + on the increasing branch (from the negatively magnetized 
state) in the final five measurement cycles is seen to be smaller than the nucleation field 
in the increasing branch of the cycles at saturation field H s . In contrast, the nucleation 
field H nt - on the decreasing branch (from the positive magnetized state) in the first five 
measurement cycles is approximately equal to the nucleation field in the decreasing branch 
of the cycles at saturated field H s . 

This asymmtery in the nucleation fields has an effect on the value of Q in equal and 
opposite bias fields, Hb and —Hb, at a given field period in sample A. In the positive bias 
field Hb, the reversal window beginning from positive saturation consists of the field intervals 
(H n _, —H + Hb) and (—H + Hb, —H p ), whereas in the negative bias field —Hb, the reversal 
window beginning from negative saturation consists of the field intervals {H n>+ , Hq — Hb) 
and (Hq — Hb,H p ). Since H nt+ < —H n _, the magnetization reversal proceeds farther in 
the bias field — Hb than in the bias field +Hb, leading to an asymmetry in the values of 
the cycle-average magnetization Q, with Q\n b > —Q\-H b - Thus, the boundary between the 
regions with negative (Q) and positive (Q) is shifted from Hb = Oe to the slightly negative 
bias field Hb ~ —1.0 Oe. 

IV. COMPARISON TO KINETIC ISING MODEL SIMULATIONS 

We now compare the experimental results to the behavior of the two-dimensional kinetic 
Ising model with nearest-neighbor ferromagnetic interactions, for which the DPT has been 
conclusively established,^ under conditions similar to those of the experiments. The energy 
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of the Ising system is given by 

E=-jY,SiS j -H{t)^2s ii (3) 

{iJ) i 

with J > the exchange constant, H(t) the time-dependent magnetic field, and Si = ±1 
the spin magnetic moment at site i. We use the Glauber transition rule, P(Si — > —Si) = 
1/(1 + e l3AE ), with the next spin to be updated chosen at random. Having verified that 
larger lattice sizes do not change the results appreciably, we use a 180 x 180 square lattice 
with periodic boundary conditions. The Curie temperature of the multilayer was estimated 
from Ref. [36| as T c « 600 K. We therefore choose the simulation temperature T = 0.5 T c , 
to correspond to the temperature T « 300 K of the experiments. As in Ref. [5|, the reversal 
dynamics are found to be MD for Hq = 0.8J at this temperature, with mean reversal time 
(r) = 59.23 ± 0.06 Monte Carlo steps per spin (MCSS). For a sawtooth waveform with 
amplitude H = 0.8J, we determine the DPT by finite-size scaling analysis 5,6 to occur at 
P c = 493 ± 2 MCSS. 

To recreate the experimental conditions as closely as possible, we simulated magnetization 
time series consisting of two saturated loops at H s = 1.6J, followed by 50 field cycles at 
Hq = 0.8J. The period values were selected from both below and above P c , from P = 368 
to 1575 MCSS. At each value of P, we collected magnetization time series for a range of bias 
fields from H b = —0.02H to 0.01H . Figures [7](a) and [7(b) show simulated magnetization 
time series, and the corresponding time series of the cycle-averaged magnetization, Qi, for 
a strong negative, weak negative, and strong positive bias field at period P = 473 MCSS, 
just below P c = 493 MCSS. As in the experimental time series in Figs. [3(a) and[3](b), the 
transition to the NESS with (Q) < in the simulated time series takes longer in the weak 
negative Hb than in the strong negative Hb , and it does not occur at all in the weak positive 
Hb- In addition, the behavior of the simulation in strong negative and weak positive Hb at 
the higher period P = 1500 MCSS, shown in Figs. [7J(c) and[7Jd), resembles the behavior of 
the experimental system at higher period shown in Figs. [3(c) and[3](d). The larger splitting 
between the average values of (Q) in Fig. E(d), relative to Fig. [7J(d), is due mainly to the 
effect of the pinning field in the experimental system, which results in an increased sensitivity 
to the bias field, as discussed later in this section. 

The similarity of the simulation data in Fig. [TJto the experimental data in Fig. [3] suggests 
that the experimental system may be near criticality just above P = 16.2 s, but this similar- 
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ity does not by itself provide conclusive evidence of a DPT. To evaluate the question more 
thoroughly, we constructed NEPDs of (Q) and a 2 (Q), measured in the NESS in simulation, 
with respect to P and H^. Two procedures were used to generate simulation data in a NESS. 
In the first procedure, we mimicked the experimental data analysis precisely, by identifying 
the part of the simulated time series which consituted the NESS using the same criterion 
described in Section IIHI In the second procedure, the simulation was initialized, for a run 
at particular values (P,Hb), such that the final saturation field before the measurement 
sequence was of the same sign as Hf,. This field-symmetric initial condition begins each 
simulation near the NESS, for both positive and negative Hb, thus eliminating the transi- 
tions seen in Fig. [3^a). We then ran 32 independent MC realizations of 50 field cycles each, 
including the final 40 cycles of each realization in the NESS data. The simulated NEPDs 
of (Q) and cr 2 {Q) using the first procedure are shown in Figs. [3(a) and[8](b), respectively, 
while those produced using the second procedure are shown in Figs. E^a) and 0(b). 

We concentrate first on the comparison of simulation data to the data from experimental 
sample A, for reasons discussed below. The simulated NEPDs of Fig. [8] are clearly similar in 
form to the experimental NEPDs in Fig. HI However, there are two significant differences, 
which can be accounted for by the presence in sample A of hard-to-reverse, residual bubble 
domains, and of a pinning field for domain-wall motion, respectively. The first difference 
is in the position of the peak of the variance a 2 (Q). In the experimental NEPD, the peak 
of u 2 {Q) is situated at (P ~ 20 s, Hb ~ —1.5 Oe), whereas in the simulated NEPD, the 
(somewhat diffuse) peak occurs at periods just above P c = 493 MCSS. While one might argue 
from Fig. [8] that the peak is centered around a slightly negative Hb, a collection of NEPDs 
(produced by different MC realizations) showed no net tendency toward either positive or 
negative Hb. As described at the end of Section IHIt the shift in the boundary between the 
regions of positive and negative (Q) values, as well as the shift of the experimental peak of 
a 2 (Q), toward negative Hb, can be explained by the asymmetry in nucleation fields during 
the measurement sequence. This asymmetry is in turn due to the presence of the residual 
bubble domains, which were magnetized positively by the final (positive) saturating field 
before the beginning of the measurement sequence. 

The second principal difference between the simulated and experimental NEPDs is the 
gradual slope of the contour lines extending from the left side of the experimental NEPD 
in Fig. H](a), as compared to the steep drop of the corresponding contours in the simulated 
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NEPD in Fig.[S^a). If one normalizes the Hb values (i.e., the y axis) by the field magnitude 
Hq in both Figs. IU(a) and[8](a), this difference in slopes becomes even more pronounced. 
Physically, the difference in slopes corresponds to a susceptibility d{Q)/dH b which falls 
off more slowly with increasing period in the experiment than in the simulation. This 
increased sensitivity to the bias field results from the pinning field in the multilayer. In the 
experiment, the entire magnetization reversal takes place within the restricted field intervals 
(— H + Hb, —H p ) and (H p , H + Hb), whereas in the simulated reversal it occurs, albeit with 
a domain wall velocity decreasing linearly with the applied field,— within the entire field 
intervals (—H + Hb, 0) and (0, H + Hb). A given bias field, expressed as a percentage of H , 
thus constitutes a larger percentage of the field interval during which magnetization reversal 
occurs in experiment than in the simulation, and so has a larger effect on (Q). 

When the second data analysis procedure was used, including averaging over 32 inde- 
pendent simulation runs, the NEPD of (Q), shown in Fig. [9](a), became symmetric with 
respect to Hb, both above and below P c . In addition, the NEPD of cr 2 (Q) in Fig.[9](b) became 
more sharply focused near Hb = at periods just above P c . This strongly suggests that the 
asymmetry with respect to Hb in Fig. H^a), and size of the peak region in Fig. H](b), would 
decrease if the second procedure described above were followed in gathering and analyzing 
experimental data. 

In comparing the simulation data to the results for experimental sample B, we note that 
in the NEPD of (Q) for sample B in Fig. [5](a), there is much less (if any) shift of the 
boundary between the region of positive (Q) and the region of negative (Q) toward negative 
bias field. Similarly, the peak of <J 2 (Q) in Figj5]^b) is just slightly further from Hb = Oe 
than the uncertainty AHb/2 = 0.5 Oe. The very minimal difference between the nucleation 
field in the increasing saturating fields and that in the increasing measurement loops in 
Fig. Mb), as compared to the sizable difference noted above for sample A in Fig. [6](a), 
suggests strongly that there were significantly fewer impurities and/or variations in film 
thickness giving rise to residual bubble domains in sample B. The NEPDs of (Q) and cr 2 (Q) 
for sample B closely resemble the part of the simulated NEPDs (Fig. [8]) above P = 500 MCSS 
P c = 493 MCSS). Thus, we tentatively place the location of the critical period at 
(P = 7.6 s, H b = 0.0 Oe), i.e., at the left edge of the NEPDs, as indicated in Fig. U This 
is also consistent with the fact that the magnetization reversal reaches m ~ for sample B 
at P = 7.6 s in Fig. EJ^b), whereas it reaches only \m\ = 0.7 for sample A at P = 8.7 s in 
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Fig. [6(a). In simulations, it has been observed 5 that the DPT occurs at a period where the 
magnetization reversal (from saturation) proceeds at least as far as m = 0. 

Since the data from sample A includes periods from both below and above our critical 
period, P c ~ 20 s, with the variance u 2 {Q) visibly decreasing in Fig. 11(b) both below and 
above this period, we have focused on sample A as our primary evidence for the observation 
of the DPT. The similarity of the experimental results from sample A to the simulated 
results, with respect to (i) the behavior of the time series Qi in various Hf,, in Figs. [3] and 
[TJ (ii) the form of the NEPDs of (Q), in Figs. IU(a) and [8(a), and (iii) most importantly, 
the existence of a well-defined peak in the NEPD of u 2 {Q) in Fig. E(b); constitute strong 
evidence for this DPT in an experimental magnetic system. The differences between the 
experimental and simulated NEPDs of (Q) and <J 2 (Q) can be consistently explained, for 
both samples A and B, in terms of the effect of hard-to- reverse bubble domains in sample 
A, as described in Section IIII| and the presence of a pinning field for domain- wall motion, 
as described earlier in this section. 

As mentioned in Section HJ the kinetic Ising model we have used does not give rise to 
a pinning field, and thus presents a simplified model of the motion of domain walls in the 
[Co/Pt]3-multilayer. Significant progress has been made in using a Preisach-Arrhenius (PA) 
model of magnetic viscosity to describe the thermally activated motion of domain walls 
over a distribution of free-energy barriers associated with variations in film thickness, grain 
boundaries, and impurities.— The PA model, while quite useful to describe phenomeno- 
logically the decay of magnetization on long time scales, does not directly incorporate the 
spatial dependence and cooperative nature of magnetization processes. In addition, recent 
work— »22i has shown that the dependence of the energy barriers (to wall motion) on the 
applied field deviates from the linear dependence^^ 1 ^ usually assumed in PA models. It 
appears likely that incorporating both spatially varying magnetization, such as that occur- 
ring in the MD reversal process (cf. Fig [2]), and a realistic description of a wide distribution 
of energy barriers will require a computationally intensive, multi-scale micromagnetic model, 
for which the DPT would first have to be established computationally. We have therefore 
chosen to model the spatially varying MD reversal process faithfully with a two-dimensional 
kinetic Ising model, and to account for magnetic viscosity effects (i.e., the pinning field) 
in our analysis and comparison to experiment. The differences between the results of the 
experiments and the kinetic Ising simulations are consistent with what one would expect 
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from the added effects of magnetic viscosity in an experimental multilayer system exhibiting 
the DPT. 

V. CONCLUSION 

We have compared the behavior of the cycle-averaged magnetization, Q, in experiments 
on a [Co/Pt] 3 multilayer, to its behavior in MC simulations of the two-dimensional kinetic 
Ising model, in which a DPT with respect to the period P of an applied alternating field has 
been well established.-^ Plots of time series of the magnetization (and of Q) in the presence 
of various bias fields Hb, as well as non-equilibrium phase diagrams of the average, (Q), 
and the variance, a 2 (Q), as functions of P and Hb, were used to characterize and compare 
the behavior of Q in experiment and simulation. The behavior was seen to be similar, with 
differences that could be clearly accounted for in terms of the known phenomena of a pinning 
field and residual bubble domains. The results, in particular the presence of a clear peak in 
the variance cr 2 (Q) in experimental sample A, provide the strongest experimental evidence 
to date for the presence of this DPT in a magnetic system. Furthermore, the results strongly 
suggest that the cycle-averaged magnetization Q serves as a dynamic order parameter in 
the [Co/Pt]3-multilayer system. 

It would be most desirable to perform further and more precise experiments on the DPT in 
perpendicular magnetic ultrathin films. Specifically, it would be useful to collect data using 
the second procedure described in Section HVl i.e., to ensure that the final saturated state 
before the measurement cycle has the same sign as the applied bias field. Also, one could 
attempt to apply a field magnitude Ho that is significantly larger than the pinning field H p , 
thereby forcing the critical period to lower values. One would expect the effect of the pinning 
field, for example in creating more gradual contours in the non-equilibrium phase diagram 
of Figs. H](a) and[Hl(a), to be reduced in this case. However, this will be experimentally very 
challenging given the hysteretic nature and response time of the electromagnet one would 
need in this case to generate the time-varying applied field sequences. 

In addition, a lower value of the critical period would enable the collection of data for 
more field cycles in each NESS. As noted above, some level of experimental fluctuations in 
the bias field is probably unavoidable. However, with a sufficiently large number of cycles 
recorded in a NESS, one could attempt to calculate values of (Q) and cr 2 (Q) using data 
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from a subset of the recorded cycles with a more restricted range of iff, values. This could 
potentially supply the precision required to investigate power-law scaling near the critical 
period. The added statistics provided by forcing the DPT to occur at a shorter period, 
as well as more precise control over the temporal stability of the bias field, could enable 
extraction of the critical exponents from the power-law behavior of cr 2 (Q) vs P — P c and Hj,. 

Finally, in an ultrathin Au/Co/Au film it was observed^ 8 that a mixture of nucleation 
at fixed ('soft', extrinsic) sites and at random (intrinsic) sites occurred. If this also occurs 
in [Co/Pt] multilayers, it could affect both the size of the fluctuations of the dynamic order 
parameter and its spatial correlations, relative to those in the kinetic Ising system, in which 
nucleation occurs only at random sites. One could examine a series of real-time images of 
the sample, such as the one shown in Fig. [21 in order to calculate spatial and time-displaced 
correlation functions to address this question. 

The results from such additional experiments could serve to further reinforce the evidence 
for the DPT presented here. To achieve an even higher degree of confidence, one would 
need to formulate an accurate, multiscale micromagnetic model of the multilayer, and to 
demonstrate the existence of the DPT in this computational model. Given the large system 
sizes needed for computational finite-size scaling analysis, such an effort will likely require 
innovations in simulation methods, as well as further increases in available computing power. 
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H (kOe) H (kOe) 

FIG. 1: Experimental data of the normalized magnetization, m, vs field H for multilayer sample 
A. Data were taken for two initial saturated loops, 49 loops at amplitude Hq = 0.366 ± 0.003 kOe, 
and one final saturated loop, (a) Period P = 16.2 s. (b) P = 38.1 s. In both cases, the bias field, 
which is defined and discussed in Section \TU\ was = +1.1 ± 0.5 Oe. 
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FIG. 2: Kerr microscope image of magnetization in a sister [Co/Pt]3 multilayer sample, taken 
shortly after the start of magnetization reversal (at m « 0.93). The area represented in the image 
had physical dimensions 0.88 x 0.62 mm 2 . 
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FIG. 3: (color online). Time series of the normalized magnetization, m, and of the cycle-averaged 
magnetization, Qi, in multilayer sample A, at two different values of the period, (a) Period P = 
16.2 s, magnetization time series, (b) P = 16.2 s, Qi vs measurement cycle index i. (c) P = 38.1 s, 
magnetization time series, (d) P = 38.1 s, Qi vs measurement cycle index %. Note that for clarity, 
only the two initial saturated cycles and the first 26 measurement cycles of the magnetization time 
series are plotted; however, the full time series of Qi (calculated from the 49 measurement cycles 
only) are shown. The bias field values corresponding to strong negative, weak negative, and weak 
positive are Hb = —3.3 ± 0.5, —0.9 ± 0.5, and +1.1 ± 0.5 Oe, respectively. 
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FIG. 4: (color online). Non-equilibrium phase diagrams (NEPDs) for sample A, showing (a) the 
average (Q) and (b) the variance cr 2 (Q) in the non-equilibrium steady state (NESS), as functions 
of the period P and bias field Hb- The portion of the magnetization time series which consituted 
the NESS was determined for each (P, Hb) according to the procedure described in Section IIII1 
Within each time series, the bias field fluctuated within a range (Hb — AHb/2, Hb + AHb/2), with 
AHb = 1.0 Oe as shown. The black dot shows the estimated location of the critical point of the 
DPT. 



0.005 



I I I M 

•C.B -.:.6 -0.4 -0,2 0.2 0.4 0,6 0.1 



CD 

o 



-0.005 




10 20 30 

Period (s) 



10 20 30 

Period (s) 



FIG. 5: (color online). Same as Fig. HJ but for sample B. 
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FIG. 6: (color online). Normalized magnetization, m, vs field, H, in (a) sample A at (P = 8.7 s, 
H b = -1.9 ± 0.5 Oe), and (b) sample B at (P = 7.6 s, H b = -0.8 ± 0.5 Oe). The upper thin 
lines (blue) show the first 5 measurement cycles, which reach (nominal) positive saturation, while 
the lower thin lines (red) show the final 5 measurement cycles, which reach (nominal) negative 
saturation. The thick black lines show two complete cycles at the saturation field, H s = 0.740 kOe. 
The nucleation fields in the measurement and saturated cycles correspond closely, except for the 
increasing branch (from negative saturation) in (a), due to the presence of positively magnetized 
residual bubble domains in sample A (see text). As the x axis range was reduced to more clearly 
show the nucleation fields, the saturation loops continue outside the field range shown. 
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FIG. 7: (color online). Time series of normalized magnetization, m, and of cycle-averaged mag- 
netization, Qi, from simulations of the L = 180 kinetic Ising model, at two different values of 
the period and in various bias fields, (a) Period P = 473 MCSS, magnetization time series, (b) 
P = 473 MCSS, Qi vs measurement cycle i. (c) P = 1500 MCSS, magnetization time series, (d) 
P = 1500 MCSS, Qi vs measurement cycle i. As in Fig.[3l for clarity, only the two initial saturated 
cycles and the first 26 measurement cycles of the magnetization time series are plotted. However, 
the full time series of Qj, calculated from the 50 measurement cycles, are shown. The bias fields 
values corresponding to strong negative, weak negative, and weak positive are —0.013 J, —0.0035 J, 
and +0.0044J. 
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FIG. 8: (color online). Non-equilibrium phase diagrams (NEPDs) for the L = 180 kinetic Ising 
model, with parameters given in Section WV\ of (a) (Q) and (b) c 2 (Q) in a non-equilibrium steady 
state (NESS), as functions of the period P and bias field H^. The data for the NESS were drawn 
from a single run of 50 field cycles, using the same analysis which was employed for the (experi- 
mental) Figs. H] and [5] and described in Section Hill The black dot shows the location of the critical 
point, (P = 493 ± 2 MCSS , Hj, = 0), as determined by finite-size scaling analysis. 
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FIG. 9: (color online). Same as Fig. El but using a different procedure to generate the data in 
the non-equilibrium steady state (NESS). Here, the data for the NESS were drawn from the final 
40 cycles and averaged over 32 independent MC simulations, each of which was initialized in a 
saturation field of the same sign as its associated H^. 
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